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Abstract 

We study the large-time behaviour of the solutions of the evolution equation in- 
volving nonlinear diffusion and gradient absorption, 

dtu- Apu+\\7u\'' = 0. 

We consider the problem posed for x <E and t > with nonnegative and compactly 
supported initial data. We take the exponent p > 2 which corresponds to slow p- 
Laplacian diffusion. The main feature of the paper is that the exponent q takes the 
critical value q = p — 1 which leads to interesting asymptoties. This is due to the fact 
that in this case both the Hamilton- Jaeobi term | Vw|* and the diffusive term ApU have 
a similar size for large times. The study performed in this paper shows that a delicate 
asymptotic equilibrium happens, so that the large-time behaviour of the solutions is 
described by a rescaled version of a suitable self-similar solution of the Hamilton- Jaeobi 
equation |VVF|p~^ = W, with logarithmic time corrections. The asymptotic rescaled 
profile is a kind of sandpile with a cusp on top, and it is independent of the space 
dimension. 
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1 Introduction and main results 



In this paper we deal with the Cauchy problem associated to the diffusion-absorption 
equation: 

dtu- ApU+lVu]" = 0, {t,x)eQ, (1.1) 
posed in Q := (0, oo) x with initial data 

u{0,x) =uo{x) >0, xGM^, (1.2) 

where the p-Laplacian operator is defined as usual by ApU := div (|Vn|^~^ Vu) . To 
be specific we take p > 2, which implies finite speed of propagation, and we consider 
nonnegative weak solutions u >0 with compactly supported initial data uq such that 

^zo G W^''°°(M^) , ^xo>0, supp (^xo) C S(0,i?o), no^O, (1.3) 

for some Rq > 0. Known properties of the equation ensure that the corresponding solutions 
will be compactly supported with respect to the space variable for every time t > 0. The 
goal of the paper is to describe in detail the asymptotic behaviour of the solutions as 
t ^ oo. 

The equation (jl.ip has been studied by various authors for different values of the param- 
eters p > 2 and g' > 1 as a model of linear or nonlinear diffusion with gradient-dependent 
absorption, see [HI [9l [TTl [121 [H] for the semilinear case p = 2, and [H [TJ [161 [2Q] for the 
quasilinear case p > 2. It has been shown that the large-time behaviour of this initial-value 
problem depends on the relative influence of the diffusion and absorption terms and leads 
to a classification into the following ranges of q: 

(i) when q > q2 := p — N/{N + 1) the large time behaviour is purely diffusive and the 
first-order absorption term disappears in the limit t — )■ oo; this is a case of asymptotic 
simplification in the sense of |21j . 

(ii) For qi := p — 1 < q < q2 there is a behaviour given by a certain balance of diffusion 
and absorption in the form of a self-similar solution, its existence being established in |2U] : 
there is no asymptotic simplification; 

(iii) for 1 < q < p—1 the two last authors have recently shown in |16j that the main term is 
the absorption term, leading to a separate-variables asymptotic behaviour, with diffusion 
playing a secondary role. We thus have asymptotic simplification, now with absorption as 
the dominating effect. 

The two critical cases q = q2 and q = qi represent limit behaviours, and as is often the case 
in such situations, they give rise to interesting dynamics due to the curious interaction of 
two effects of similar power. Such situations usually lead to phenomena called resonances in 
mechanics, with interesting non-trivial mathematical analysis. Such interesting behaviour 
has been shown in particular in [11] for q = q2, in the linear case p = 2, with the result 
that logarithmic factors modify the purely diffusive behaviour found for q > q2- A similar 
situation is expected to be met when p > 2 and q = q2- 

We devote this paper to describe the other limit case q = qi = p — 1 when p > 2, the 
latter condition guaranteeing that q > 1. In that case the diffusion and the first order term 
have similar asymptotic size and logarithmic corrections appear in the asymptotic rates. 
The mathematical analysis that we perform below is strongly tied to a good knowledge of 
the expansion of the support of the solution, or in other words, the location of the free 
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boundary, which happens to be approximately a sphere of radius \x\ Clogt for large 
times t. From now on, we assume that 

q = qi=p-l. 

1.1 Bounds in suitable norms 

Studying the large time behaviour of solutions and interfaces of our problem relies on 
suitable and very precise estimates. The time expansion of the support and the time 
decay of solutions to the Cauchy problem ()l.ip - (jl.2p . with non-negative and compactly 
supported initial data have been recently investigated in [7j. The following results are 
proved: 

Proposition 1.1. Under the above assumptions on the equation and data, the Cauchy 
problem has a unique non-negative viscosity solution 

u G SC([0, oo) X M^) n L°°(0, oo; W^'^{mJ^)) 

which satisfies: 

< u(t,x) < lltiolloo, (t,x)EQ, (1.4) 

||Vn(t)||oo < llVnolloo, t>0, (1.5) 

supp {u{t)) C -8(0, Ci logt) for all t>2, (1.6) 
together with the following norm estimates 

\Ht)\\i < (logt)(P(^+i)-2^-^)/(P-2) j^^^ii ^>2, (1.7) 

\\u{t)\\oo < C2 t-^/^^'-^) (logt)(P-^)/(P-2) j^^^ii ^^2, (1.8) 

||V7/(t)||oo < C2 (logt)i/(P"2) for all t>2, (1.9) 

for some positive constants Ci and C2 depending only on p, N , and uq. 

Here and below, BC{[0, 00) x M^) denotes the space of bounded and continuous functions 
on [0,00) X and || • \\r denotes the L'"(]R^)-norm for r £ [l,oo]. As we shall see, 
these bounds will be very useful in the sequel. The well-posedness of (|l.l|) - (|1.2p and the 
properties (|1.4p . (|1.6p . and (|1.7p are established in [71 Theorems 1.1 & 1.6, Corollary 1.7], 
while (|1.8p and (|1.9p follow from (|1.7p and [71 Proposition 1.4]. We will also use the 
notation r+ = max{r, 0} for the positive part of the real number r. 

1.2 Main results 

We describe next the main contribution of this paper. As already mentioned, our goal is to 
study the asymptotic behaviour of the solution u of the resonant problem (jl.ip with p > 2 
and q = p — 1, and with compactly supported and nonnegative initial data. Moreover, 
since the equation has the property of finite speed of propagation, it is natural to raise the 
question about how the interface and the support of the solution expand in time. We also 
answer this question in the present paper. 

Asymptotic behaviour. The main result is the following: 
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Theorem 1.1. Let u he the solution of the Cauchy problem ()l.ip - ()1.2p . with uq as in 
(|1.3p . Then, u decays in time like 0{t~^/^'P~'^\logt)^P~^^/^'P~'^^) and the support spreads 
in space like O(logt) as t ^ oo. More precisely, we have the limit: 



lim sup 



0, (1.10) 



(logi)(P-i)/(P-2) ' ' ' V log t 

OTi/i precise constant 

Cp = (p_2)V(p-2)(p_l)(p-l)/(p-2), 

In the proof, the expression of the asymptotic profile is obtained after a comphcated time 
scahng of u and x in the form of uniform hmit 

fl/(p-2) 

n(t,x) ^ (p-2)-P/(P-2) P^((p-2)2;/logt), (1.11) 



(logt)(P-i)/(P-2) 
where the function 

/„_2 \(p-i)/(p-2) 

is the unique viscosity solution to the stationary form of the rescaled problem, which is: 
\VW\P-^-W = in 5(0,1), H^ = on dB{0,l), > in 5(0,1). (1.13) 

Let us notice that, as usual in resonance cases, the limit profile is not a self-similar solution, 
but it introduces logarithmic corrections to a self-similar, separate- variables profile (which 
in our case is t~^^^~'^\p — 2)~'P^^^~'^^W{{p — 2)x)). The uniqueness of W as viscosity 
solution of (jl.l3p is very important in the proof and follows from |13j . 

In consonance with (jl.lOp . we show that the shape of the support of u{t) gets closer to a 
ball while expanding as time goes by. This is in sharp contrast with the situation described 
in [16] for (jl.ip in the intermediate range q € {l,p — 1), p > 2 where the positivity set 
stays bounded and can have a very general shape. When q = p — 1, the diffusion thus 
acts in three directions: the scaling is different, the support grows unboundedly with time, 
and the geometry of the positivity set simplifies. Another remarkable consequence of the 
interplay diffusion-absorption is the fact that the asymptotic profile is radially symmetric 
and does not depend on the space dimension. 

We devote Section U] to the proof of Theorem ll.il For the proof, we use a precise estimate 
for the propagation of the positivity set, that is described below. Another tool is the 
existence of a large family of subsolutions having a special, explicit form and allowing for 
a theoretical argument with viscosity solutions to finish the proof. 

Propagation of the positivity set. We denote the positivity set and its maximal 
expansion radius by 

Vuit):={xeR^ : u{t,x)>0}, 7(t) = sup{|x| : x G P„(t)} (1.14) 

respectively. Then: 

Theorem 1.2. Under the running notations and assumptions, we have: 

lim ^ = (1.15) 
t^oo logt p-2 

Moreover, the free boundary of u has the same speed of expansion in any given direction 
w G with \uj\ = 1. 
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In fact, we give more precise estimates for the expansion of the positivity region, ob- 
tained via comparison with some well-chosen traveling waves. The proof of Theorem 11.21 
is performed in Section [3l 

Two scalings. In order to prove the two theorems, we have to perform two different scaling 
steps. The first scaling, described in formula ()2.2p below, is the natural one corresponding 
to standard scaling invariance; such a scaling has also been used in [16] in the case q G 
(l,p — 1) to obtain the correct scale of the solutions. But for g = p — 1, we observe that 
a phenomenon of grow-up appears, which is typical for resonance cases: the effect of the 
resonance implies that the rescaled solution does not stabilize in time; on the contrary, it 
grows and becomes unbounded in infinite time. That is why we need a second scaling, 
given by the new functions w and y defined in (j4.ip and (j4.2|) . which is less natural but 
turns out to be adapted to our problem: it takes into account the logarithmic corrections 
(suggested by the a priori estimates of Proposition ll.H which turn out to be sharp), and 
it is adapted to the size of the grow-up phenomenon; thus, in the rescaled variables we can 
describe the real form and behaviour of the solution. 



2 Scaling variables I 

We recall that p > 2 and q = p — 1. We introduce a first set of self-similar variables; we 
keep the space variable x and introduce logarithmic time 

r:=-^ln(l + (p-2)t), (2.1) 

as well as the new unknown function v = v{t, x) defined by 

u{t,x) = {l + {p-2)ty^/^P~^K{T,x) , (i,x) G [0,oo) xM^. (2.2) 
Clearly, v solves the rescaled equation 

drV - ApV + \Vv\'' - V = , {T,x)eQ, (2.3) 
with the same initial condition 

v{0) = uo, x€R^ . (2.4) 

We next translate the a priori bounds (jl.7p . (jl.Sp . and (jl.9p in terms of the rescaled 
function v: there is C3 > depending only on p, N, and uq such that 

+ ^T77;;z^ ^ '^3 for r>l. (2.5) 



T-(p(7V+l)-27V-l)/(p-2) ' ^(p-l)/{p-2) ' ^l/(p-2) 

2.1 The positivity set: time monotonicity 

We define the positivity set Vv{t) of the function v at time r > by 

P^(r) := {x G : v{t,x) > . (2.6) 
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Proposition 2.1. For ti € [0, oo) and T2 G (ti,oo) we have 

V,{ti) <^V,{t2) and lJp,(r)=M^. (2.7) 



T>0 



/n addition, for each x G M there are Tx > and e^^ > such that 

v{t,x) > Ex t(p-^)/(P'2) T>Tx. (2.8) 

The proof relies on the availabihty of suitable subsolutions which we describe next. 
Lemma 2.1. Define two positive real numbers Rp and Tp by 

Rp--=^^^^ and Tp:=^-^^ {2 + 2^-HN+p-2)) . 
If R £ (0, Rp] and T > Tp, the function sr^t given by 

is a (viscosity) subsolution to ()2.3p . 

Proof. We have sr^t{t,x) = (T + Tfv-~i)/{p~2) ^^^^ ^-^.j^ ^ _ ^/(^j. _^ g^^^^ ^(^^^ _ 
(p - 2) (ii^ _ |^|2^(p-i)/{p-2) /(^^^ _ gij^(,g j3-l>p-2>0, we observe that a and 
iVdl^'^^Vcr both belong to ^^(IR^). Therefore, 

L(t, x) :=R{T + rr^P~mP~2) _ ^^^^^^ ^ | Vs«,Tr-' - ^/j.t} 

is well-defined for (r, x) € [0, 00) x and 

l(t,x) = ^ |^a(o-e-va(o-A,c7(e)} + i2 iv<7(e)r^-i?<7(e) 

- (^2_,^,.)(.-i)/(.-2) (l + 2-H^ + p-2)f^)} 

+ (^2 _, {2-^f^-^-^} 



< {R^-mt'^^'^-'' j^^±i:!iN±P^^,p^,^_P_2 

^ /r2 |,|2Up-1)/(P-2) / 2 A 2P-l(p-l) ler^X ^ 



We next note that 



p - 2 - i^-i - 2 ' 
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so that the last term of the right-hand side of the previous inequahty is bounded from 
above by 2 {R"^ - |^|2)(p-i)/(p-2) ^^3^^^ Consequently, owing to the choice of R and T, 

L(r..) < (R^ - Itf )r"^'--' I ^ + ^'-(jy + P - ^) , - ^ . J. 

< 0, 
whence the claim. □ 

Proof of Proposition 12.11 (i) Fix ri > and xi € Vv{ti). Owing to the continuity of 
X I — > v{ti,x) there are 5 > and ri > such that v{ti,x) > 6 ior x ^ B{xi,ri). Take 
now R > small enough such that R < minjri, Rp} and satisfying 

R < and ^ (Tp + ri)^P-mP-2) rp/{p-2) < ^ 

Tp + Tl p-1 

the parameters Rp and Tp being defined in Lemma |2. II Then we have SR^TpiTi,x — xi) = 
0<v{ti,x) if \x — xi\ > R {Tp + Ti), while 

SR,TAn,^- ^1) < J^^^ {Tp + r,)^P~mP^^) R(2p-mp-2) < ^ < 

if |a; — xi I < i? (Tp + ti) as -R(Tp + ri) < ri. Moreover, if r2 > ri, r € [ti,T2] and x G 
dB{xi, R{Tp+T2)), then sj^^Tp(''"5 2; — xi) = < v{t, x). Recalling that s_R,Tp is a subsolution 
to (j2.3p by Lemma [27n we infer from the comparison principle that SR^Tp {t, x—xi) < v{t, x) 
for {t,x) e [ti,T2] X B{xi,R{Tp + T2)). As s_r,Tp(t, x-xi) = < i)(r,3;) for r G [ri,r2] and 
X ^ -B(xi, i?(Tp + T2)) we actually have sr^Tp{t, x — xi) < v{t, x) for (r, x) € [n, T2] x M^. 
Since T2 > ti is arbitrary and neither R nor Tp depend on T2, we end up with 

SR,Tp{T,x - xi) <v{t,x) , (r,x) e [ti,cx)) X . (2.9) 

A first consequence of (|2.9|) is that, if T2 > ti, then i;(T2,xi) > sr^Tp{t2,0) > so that xi 
also belongs to Vv{t2)- 

Next, given x € M^, we have x € B{xi, R{Tp + r)) for r large enough and it follows from 
(|2.9p that v{t,x) > sr^Tp{t,x — xi) > for r large enough. Consequently, x belongs to 
Vv{t) for r large enough which proves the second assertion of ()2.7p . 

(ii) Consider xq € M^. According to ()2.7p there is tq large enough such that xq € Vv{tq). 
Arguing as in the proof of (|2.7p . we may find tq small enough (depending on xq) such that 
Sro,Tp{T,x — Xq) < v{t,x) for (t, x) G [to,oo) X M.^ . Consequently, 

y(r rnl > ^"^ (T +r^(P-l)/(P-2) (2p-2)/{p-2) p/(p-2) (p_i)/(p_2) 



which gives the lower bound ()2.8p . □ 

Corollary 2.1. Assume that uq{0) > 0. T/ien i/iere is r* > suc/i i/iai 

o^ / 1^12 \ (p-i)/(p-2) 

l^'-^i n , _\(p-l)/{p-2) /^„2 m 



Proof. Arguing as in the proof of (j2.7p and using the positivity of uq{0), we may find 
r* > small enough such that Sr,^TpiT,x) < v{t,x) for (r, x) G Q. Since Tp > 1, the 
previous inequality implies (j2.10p . □ 
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2.2 Eventual radial symmetry 

We prove the following classical monotonicity lemma, see [31 Proposition 2.1] for instance. 
Lemma 2.2. If x G and r > satisfy \x\ > 2Rq and r < |x| — 2Rq. Then, 

v{t,x) < inf v{T,y) for t>0. (2-11) 

\y\=r 



Here, Rq is radius of the initial ball defined in (jl.3p . 

Proof. The proof relies on Alexandrov's reflection principle. Let (x, r) G x (0, oo) 
fulfil the assumptions of Lemma 12.21 and consider y € such that \y\ = r. Let H be the 
hyperplane of points of which are equidistant from x and y, namely 



H:=->zeR^ : (z-^^,x-y) =0 



Litroducing 



H_:=izeR^ : (z-^^,x-y) <0 



and 



v{t, z) := V { t, z - 2 (z- x - y^ ) , {t,z) eQ , 



it readily follows from the rotational and translational invariance of ()2.3p that v also solves 
dOjl . In addition, y € and P„(0) C B{0,Ro) C by ([L3]). Now, on the one hand, 
if z € then 

x+y \ x-y 



2 / |x — 

and v(0, = < f (0, z). On the other hand, \i z ^ H = dH^ and r > 0, we clearly have 
v{t,z) = v{t,z). We are then in a position to apply the comparison principle to (|2.3|) on 
(0, oo) X and conclude that 

v{t, z) < v{t, z) , (r, z) G [0, oo) X . (2.12) 

Recalling that y G H^, we infer from (|2.12|) that v{t, y) > v{t, y) = v{t, x) for r > which 
is the expected result. □ 

Remark 2.1. Although Lemma \2.S\ will not he used in the main proofs, this is an inter- 
esting result for the qualitative theory, since it shows that the dynamics symmetrizes the 
solution. 



3 Propagation of the positivity set 

We next turn to the speed of expansion of the positivity set Vv of v and put 

^(r) := sup{|x| : xGPt,(r)}, (3.1) 

so that Vv{t) Q B{0,g{T)) for r > 0. The purpose of this section is to prove that the 
expansion speed g{T) of Vv{t) is asymptotically equal to r, in other words, 

hm ^ = 1, 

T— >oo T 
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and, more precisely, to prove Theorem 11.21 

The proof rehes on the existence of "nice" travehng wave solutions of (|2.3p , which may 
be used as subsolutions and supersolutions for the Cauchy problem ()2.3p - ()2.4p . The con- 
struction of such traveling waves is inspired on the technique used in the so-called KPP 
problems, [H], which has developed a wide literature; see e.g., [2], [22] for applications to 
porous media, and [18] for blow-up problems. We thus begin with a phase-plane analysis, 
proving the existence of the desired traveling waves. 



3.1 Traveling wave analysis for = 1 

We look for traveling waves of the form 

v{t, x) = f{z), Z = X — CT, c > 0, 



solving (j2.3p in dimension = 1. Then, the profile / solves the ordinary differential 
equation: 

-c/'-(i/'rY)'+i/'r'-/ = o. (3.2) 

We are actually only interested in traveling waves which present an interface, that is, / 
vanishes for z sufficiently large. As we shall see below, the profile / is non-monotone in 
general, but is nonnegative and decreasing near the interface. We transform (j3.2p into a 
first order system, by introducing the notation U = f and V = —f. We arrive at the 
following system 

{p - i)\vr^u' = -{p - i)|yrv, 

{p - i)\v\p-^v' = -cv - \v\p-^ + u, ^ ' 

where, for the orbits, the term (p — 1)|1/|*'~^ in the right-hand side has no influence (since 
we work with dV/dU) and can be ignored after a change of the time variable. We perform 
next the phase-plane analysis of the system (j3.3p . 



Local analysis in the plane. The system ()3.3p has a unique critical point, P = (0,0), 
and the Jacobian matrix J(0, 0) at this point is given by 



with eigenvalues Ai = and A2 = — c, and corresponding eigenvectors are ei = (c, 1) and 
62 = (0, 1). By a careful analysis, we notice that the center manifold in P is tangent to ei, 
and is asymptotically stable. It follows that P is a stable node for every c > 0. There is a 
unique orbit entering P and tangent to 62, forming the stable manifold; its local behaviour 
is U{z) ~ C(-z)(P"^)/(P"2) as z 0, hence this orbit contains all the traveling waves with 
velocity c and having an interface. By standard theory (see, e.g., [17]), all the other orbits 
approach the center manifold, tangent to ei, and present an exponential decay, but no 
interface: U{z) ~ e~'^^ as z — > 00. 

Local analysis at infinity. We investigate the behaviour of the system when U is very 
large. For monotone traveling waves, we make the following inversion of the plane: 

.4, w-J^. 
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and we are interested in the local behaviour near Z = 0. After straightforward calculations, 
p.3p becomes the new system: 



( Z' = z(2p-3)/{p-i)vF|VF|~(P~2)/(p-i)^ 

{ W' = z(P-2)/(p-i)|t^|p/{p-i) - cZ(P-2)/(P-i)VF|VF|-(P-2)/(P-i) + 1 - \W\. ^^'^^ 

We find two critical points with Z = 0, namely Qi = (0,1) and Q2 = (0,-1). We wiU 
analyze only Qi, i.e. the decreasing traveling waves. Let us also remark that, in the second 
equation of ()3.4p . the terms with Z are dominated by 1 — near Qi and Q2, hence we 
can study the local behaviour by using the approximate equation only with 1 — \W\ in 
the right-hand side. The linearization near Qi has eigenvalues Ai = and A2 = —1, and 
the center manifold, which is tangent to the line = 1, is unstable. Hence, the point 
Qi behaves like a saddle, and the orbits which are interesting for our study are the orbits 
going out of Qi. These orbits are tangent toW = l, and in the original system they satisfy 
U r-^ V^~^, hence, by integration, 

C/(z) ~ asz^-00, 

and are decreasing. The local analysis around Q2 is similar, but not interesting for our 
goals. 

Let us notice that not all solutions passing through a point in the first quadrant come 
from Qi. Indeed, the orbits touching the curve U = cV + V^~^ change monotonicity as 
functions V = V{U), hence they have previously reached the axis ^ = 0, meaning a change 
of monotonicity as / = f{z), and they enter through this change in the first quadrant. 
Analyzing the curve U = cV + V^~^, we observe that it connects in the phase-plane the 
points P = (0,0) and Qi, being tangent in Qi to the axis W = 1. In particular, there 
exist non-monotone solutions, and this is the object we are interested in. 

Global behaviour. This is now not difficult to establish, by merging the previous local 
analysis with the following important remarks: 

(a) The evolution of the system ()3.3p with respect to the parameter c is monotone. Indeed, 
we calculate: 

d^ [dU ) ~ (p- 1)|T/|P"2 ^ 

(b) There exists an explicit family of traveling wave solutions with speed c = 1: 

^■^(^)=(^) (K - z)t'^/^^"'^ ■ K>0, (3.5) 

This function is obviously decreasing and presents an interface at z = K. It is immediate 
to check that this orbit satisfies U = V^~^, hence it comes from the point Qi along the 
center manifold of it, and it enters P, being the unique orbit entering P and tangent to 
the eigenvector 62 = (0, 1) (unique for c = 1), as discussed above. 

(c) Moreover, the vectors of the direction field of ()3.3p over the curve U = V^~^ (which 
gives the explicit orbit (|3.5p ) have the same direction. Indeed, the normal vector to this 
curve is (1, —{p — 1)VP~'^) and we calculate: 

(1, -{p - l)yP-2) . {-{p - l)VP-\ -cV - VP-^ + U) = {p- l)(c - i)vp-^. 
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For c = 1 we obtain the explicit trajectory, and for c < 1, the above scalar product is 
negative, hence all these vectors have the same direction, contrary to (1,— — 1)^^"^). 
For c > 1, all these vectors have the same direction as V. 

Since we are interested only in traveling waves with an interface, we analyze only the 
unique (for c fixed) orbit entering P = (0, 0) tangent to 62 = (0, 1). For c = 1, it is explicit 
and connects P and Qi in the first quadrant. We draw the phase-plane for c = 1 in Figured] 
below; it is clear that the explicit connection will not change sign and monotonicity. 




Figure 1: Phase portrait around the origin for c = 1. Experiment for p = 3, N = 2. 

By remarks (a) and (c) above, it follows that for c < 1, this unique orbit disconnects from 
Qi, hence it should cross at some point the curve U = cV + VP~^ (which still connects 
P = (0, 0) and Qi); as explained before, this orbit previously had a change of sign (crossing 
the axis U = 0) and then a change of monotonicity (crossing the axis V = 0). In particular, 
we can say that the explicit orbit ()3.5p is a separatrix between the monotone and the non- 
monotone orbits. We draw the local phase portrait for c < 1, around the origin, in Figure[2] 
below. We gather the discussion above in the following result. 

Lemma 3.1. (i) For any c € (0, 1) and K > 0, there exists a unique traveling wave solution 
f c,k{^) — f c,k{^ ~ ct) of (j2.3p in dimension N = 1, having an interface at z = K (that 
is, fcxi^) — fo''^ z > K) and moving with speed c. In addition, fcxi^) — /co(^ ~ ^) 
for z £R. 

(a) For c = 1 and for any K > 0, there exists a unique nonnegative traveling wave 
fi,K{z) = fi,K{x — t) of ()2.3|) in dimension N = 1 with interface at z = K, having the 
explicit formula: 

/i,i.(x-r)=(^) (/^ + r-x)r)/(-^). (3.6) 

Here again, fi^xiz) = fi,o{z — K) for z G M. 

(Hi) For any c > 1 and K > Q, there exists a unique traveling wave solution fc^K = 
fc,K{x — ct) of ()2.3p in dimension N = 1 with interface at z = K and moving with speed 
c. Moreover, fc^K is nonnegative and non-increasing, and fc,K{z) = fcfi{z — K) for z S M. 
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-3 r 

-2 







Figure 2: Phase portrait around the origin for c < 1. Experiment for p = 3, = 2, 
c = 0.9. 



Compactly supported subsolutions for < c < 1. We are looking for nonnegative and 
compactly supported subsolutions traveling with any speed < c < 1. These subsolutions 
are constructed in the following way: from the analysis above, we know that, given c € 
(0, 1) and K > 0, there are two points Zc,k S (— and Zc^k £ izc,K,K) such that 

Zc,i^ := inf {z G (-00, : fc,K>^ in iz,K)}>-oo, 

and 

fc,K>0 in {zc,k,Zc,k) and fc,K<^ in {zc,K,K). 
We then define _ 

f f^\-j fc,Kiz), for Zc,K <z<K, 

^^'""^'^-[o, elsewhere. ^^'^^ 

In other words, we consider the positive hump of the graph of fcj< located between its 
last change of sign and the interface. It is immediate to check that fc,K is a compactly 
supported subsolution to ()2.3p in dimension = 1, and that it has an increasing part until 
reaching its maximum at Zc,k, and then decreases to the interface point K. The notation 
fc,K will designate in the sequel these subsolutions if < c < 1 and the solutions to ()2.3p 
in dimension = 1 given by Lemma l3. II if c > 1. 



3.2 Construction of subsolutions in dimension > 1 



We turn to equation ()2.3p posed in dimension A^ > 1 for which we aim at constructing 
some special subsolutions having an interface that moves out in all directions with a given 
velocity c < 1. The construction is based on the traveling waves /c,x identified in the 
previous subsection. The first attempt is to try the form V{t, x) = fc^Ki\x\—CT), c G (0, 1), 
which satisfies: 

drV - ApV + \VV\P-^ - V 



< hf> \v-'if> 

\Za^\jc,k\ Jc,K ■ 
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Thus, V is a subsolution of ()2.3p in the region of Q where > 0. We therefore have to 
modify the profile in the decreasing part of fc.K and we proceed as follows. 

Traveling wave solutions to a modified equation in dimension = 1. For a G 

(0, 1/2), we consider the following perturbation of ()2.3p : 

5.c-54ia,cr'5,c) + i5xcr'-a|5xcr'5,c-c = o, {t,x) g (o,oo) xm, (3.8) 

and look for traveling wave solutions C(t, x) = g[x — ct). Then, g solves 

- eg' - {wr'g'y + wr' - « wr'g' -9 = 0. m 



The phase-plane analysis for (|3.9p is similar to that of (|3.2p . with the difference that an 
extra term —a \ V\p~'^V appears in the right-hand side of the second equation in (j3.3p . This 
is only reflected in the analysis at inflnity, where the point Qi changes into (0, 1/(1-1- a)) 
and the explicit separatrix is obtained for c = 1/(1 -|- a) < 1. In particular, we have the 
following analogue of Lemma 13.11 (i) . 

Lemma 3.2. For any a > sufficiently small, c (z (0, 1/(1 -|- a)) and K >0, there exists 
a unique traveling wave solution gc,K,a{z) = gc,K,aix — ct) of ()3.8p having an interface at 
z = K and moving with speed c. In addition, gc,K,a{z) = gc,o,a{z — K) for z G M and there 
are two points Zc^K,a G ("OO, ET) and Zc^K,a £ izc,K,a, K) such that 

Zc,K,a ■■= inf {z G (-oo, K) : gc,K,a > in {z, K)} > -oo , 

and 

9c,K,a>0 in iZc,K,a,Zc,K,a) and g'c,K,cx<^ {Zc,K,a,K). 

Setting 



Mc,a := sup {5c,0,a(^)} , 



we notice that 



Zc,K,a = Zc,Q,a + K , Zc,K,a = Zcfi,a + K , SUp {gc,K,a{^)} = Mc,a ■ (3.10) 

ze[zc,K,c,,K] 

If we put now V{t,x) = gc,K,ai\x\ — ct), we calculate and find that 

drV - A,V + \VVr' - y = (a - ^) (|9U,«r'9U«) (1^1 - ' 

and it is a subsolution where g'^j^^<0 and a > (A^— l)/|x|. Matching these two conditions 
turns out to be possible as we show now. 
Fix c G (1/2, 1) and ac := (1 — c)/(l -|- c) and define 

f 2(iV-l) ge,0,a. K - 1) r-^m 

ro(c) := max <^ 2zc,o,a,, > > , (3.11) 

[ ac c } Uc 

the point ic,o,«c ^ (— co,0) being defined in Lemma [3.21 Then c < 1/(1 -|- Uc) and, for 
K > 0, T > ro(c), and \x\ > Zc,K,ac + ct = Zcfi,ac + K + CT, we have 

iV-1 N-l 2{N-l) 

\A Zcfi^a^ -h Cro(c) 22:c,0,ac + '^o(c) 
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and 



9c,K,aS\^\ - ct) < if Zc,K,ac+ CT < \x\ < K + CT , 
9'c,K,aS\^\ ~ ct) = if K + CT<\x\. 

Consequently, for c € (1/2, 1), Oc = (1 — c)/(l + c), and A' > 0, the function V defined by 
,x) — gc,K,ac{\^\ — cr) is a subsolution to ()2.3p for r > to(c) and |x| > ^c,/c,qc + ct. 
Observing that any positive constant is a subsolution to (|2.3|) . we construct a compactly 
supported subsolution Vc^k to ()2.3p by setting 

{Mc,ac if < |x| < Zc,K,ac + CT , 

(3.12) 
gc,K,aA\A - ct) if |x| > Zc,K,a, + CT , 

for r > To(c). It is easy to check that the function fc,x is a subsolution to (|2.3p in 
[ro(c),oo) X M^. It will be used for comparison from below, as indicated in the next 
subsection. 



3.3 Proof of Theorem [IH 

We conclude the proof of Theorem 11.21 bv a comparison argument, using the subsolutions 
and supersolutions constructed in the previous subsections. Before that, we identify a class 
of solutions of ()2.3p that is representative for the general solutions. 

We say that a function V = V{t,x) is radially non-increasing if V{t, ■) is radially sym- 
metric for all r, and it is non-increasing in the radial variable r := For example, 
the subsolutions Vc^k are radially non-increasing. The next results show that the class of 
radially non-increasing solutions of ()2.3p is sufficient for our aims. 

Lemma 3.3. Let uq = UQ{r) be a radially non-increasing function satisfying (|1.3p . Then, 
the solution v of (|2.3p with initial condition uq is also radially non-increasing. 

Proof. The radial symmetry of the solution v follows from the invariance of the equation 
(j2.3p with respect to rotations. We write now the equation satisfied by ^ = drV, obtained 
by differentiating (|2.3p with respect to r: 

- dhi^r^) - ^^driicr^o + ^^icr'e + {p- mr'^Or^ - e = o, 

which is a parabolic equation (of porous medium type) and satisfies a maximum principle. 
Since is a solution to the above equation, the derivative ^ = drV remains nonpositive if 
it is initially nonpositive and it follows that v is radially non-increasing. □ 

We are now in position to end the proof of Theorem ll.2l for radially non- increasing initial 
data. More precisely, we have the following upper and lower bounds for the edge g (r) 
defined in ()3.ip of the support of v{t). 

Lemma 3.4. Let uq = UQ{r) be a radially non-increasing function satisfying ()1.3p and 
denote the solution of (j2.3p with initial condition uq by v. For any c G (1/2,1), there 
exists Ti(c) > such that, for any r > ti(c), we have: 

1 + c{t-ti{c))<q{t)<Ro + ^ ||no||(P-2)/(P-i)+r. (3.13) 

p — 2 

In particular, we obtain that q{t)/t 1 as r —)• oo. 
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Proof. The upper bound follows by comparison with the explicit traveling wave solutions 
p.6p . More precisely, we define 

Ri:=Ro + ^ \\no\t-'^^^'-'^ (3.14) 
p — 2 

and consider the function v{t, x) = (xi — t), which is a solution of (|2.3|) by Lemma [3. 1[ 
If X = (xj)i<j<Ar € is such that xi > Rq, then |x| > Rq and uo{x) = < v{0,x) while, 
if xi < Rq, 

/„_ 9\ {P-l)/(P-2) , , 

no{x) < llnolloo < (^^j (R^-Ro)^-'^/^^-'^ 

/r, 9\ (P-l)/(P-2) , ^ 

< (^) (i?i-xi)(^-^)/(^-^)=iJ(o,x). 

The comparison principle then entails that v{t,x) < v{t,x) for (r, x) S [0, oo) x M^, 
from which we conclude that P^(t) C {x G : xi<Ri+t}.0 wine to the rotational 
invariance of (j2.3p . we actually have Vi,{t) C |a; G : {x,uj) < -Ri + r} for every cj G 
and T > 0, and thus 

V^{t)CB{0,Ri + t). (3.15) 

The lower bound follows from comparison with the subsolutions constructed in ()3.12p . 
Fix c G (1/2,1) and put ri := 1 + cro(c), To(c) being defined by ()3.1ip . Since v{t) is 
radially non-increasing for all r > by Lemma [3.31 we infer from Proposition 12. II that, for 

X G B{0,ri) and r > T^^, 

Define ti(c) by 



ri(c) := max < To(c),r,.,, 



^c,{l-c)/(l+c) A 



SO that the previous inequality and the properties of t)c,i defined in ()3.12p guarantee that 

v{ti{c),x) > Mc,(i_c)/(i+c) > Vc,iiTo{c),x) , X G B{0,ri) . 

Since Vc^i{tq{c), x) = for x i?(0, ri), we also have f(Ti(c),x) > t>c,i(ro(c), x) for x ^ 
B{0,ri). Recalling that Vc^i is a subsolution to ()2.3p in (ro(c),oo) X R^, we infer from the 
comparison principle that 

v{t + Ti{c),x) > Vc,iiT + To{c),x), {T,x)eQ. (3.16) 

Consequently, v{t + ti(c), x) > if x G i?(0, ri + cr), whence 

S(0,1 + c(t + to(c)-ti(c))) cP„(t), r>ri(c). (3.17) 

This readily implies that 

q{t) > 1 + c(t + To(c) - Ti(c)) > 1 + c(r - ri(c)) , r > ri(c) . 
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In particular, we deduce from ()3.15p and ()3.17p that 

liminf ^ ^ ^ > c for any c G (1/2, 1) and limsup ^ ^ ^ < 1 , 

r^oo r T^oo T 

which impUes that q{t)/t 1 as t ^ oo. □ 

Rephrasing the rescahng and coming back to the notation with t = {e^P~'^^'^ — l)/(p — 2) 
and j{t) = q{t), we find the result of Theorem 11.21 for radially non-increasing inital data. 
The extension to arbitrary initial data satisfying (jl.3p is performed in Section [5j Moreover, 
we notice that the speed is the same in any direction u € as stated. 

4 Proof of Theorem 11.11 
4.1 Scaling variables II 

According to Proposition 12. H as r — )• oo the solution v to ()2.3p . ()2.4p expands in space and 
grows unboundedly in time. In order to take into account such phenomena, we introduce 
next a further scaling of the space variable 

2/:=T^, (4.1) 
1 + r 

together with the new unknown function w = w{t, y) defined by 



v{t, x) = (1 + TfP-'^I^P-'^ w (^r, ] , (r, x) e [0, oo) x M^" . (4.2) 

It follows from (12.31) and (12.41) that w solves 



drW---^(ApW + yVw-- — Iw) +\Vw\P~^ -w = , {T,y)€Q, (4.3) 
l + ry p — I J 

with the same initial condition 

w{<d) =uo, y G . (4.4) 

Throughout this section we assume that uq is radially non-increasing besides (|1.3p . In 
particular, no(0) > 0. We gather several properties of w in the next lemma. 

Lemma 4.1. There is a positive constant C4, depending only on p, N, and uq such that 

||w;(T)||i + ||w;(r)||oo + ||Vu;(T)||oo <C4, r>0, (4.5) 

w{r,y) > ^ {rl - \y\'il-'^ ' ^"-'^ , {r,y) e Q , (4.6) 
the radius r^, being defined in Corollaru \2.1\ Moreover, 

P^(t) := {y € : u;(t, y) > O} C i? (^0, 1 + (4.7) 

for T > where Ri is defined by (|3.14p . In addition, for any c € (1/2, 1), we have 

B ( 0,c- ) C V^{t) for r > ti(c) , (4.8) 



1 + ^, 

the time ti(c) > being defined in Lemma\3.4 
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Proof. The estimates and (|T6]l readily follow from ([23]) and ([2J0]) . while (fO]) is 
a consequence of (|3.15p . The assertion about the ball .6(0, c — ri(c)/(l + r)) follows from 

(IMZD. □ 



At this point, (|4.3p indicates that ii;(r) behaves as r — )• oo as the solution w to the 
Hamilton- Jacobi equation drW + iV'iZ;!^"^ — zZ; = in Q which is known to converge to a 
stationary solution uniquely determined by the limit of the support of w(t) as r — > oo, 
see, e.g., [15] Theorem A. 2]. As an intermediate step, we thus have to identify the limit 
of the support of t(;(r) as r — oo. Thanks to (|4.7p . we already know that it is included 
in B{0, 1) but the information in (|4.8p are yet too weak to exclude the vanishing of w{t) 
outside a ball of radius smaller than one. To complete the proof of Theorem 1 1.1 1 for radially 
non-increasing initial data, we show first that the asymptotic limit is supported exactly in 
the ball -6(0, 1). Then we use a viscosity technique, the same that has been used in the 
previous paper [16] to establish the convergence to the expected stationary solution. 



4.2 Proof of Theorem [Hit = 1 

We first consider the one-dimensional case = 1 and divide the proof into several technical 
steps. 

Step 1. A special family of subsolutions. Given c G (1/2,1), we have 

v{t,x) > Vc,i{t + To(c) - ri(c),x) , {t,x) G [ti(c),oo) X M, 

by ()3.16p . the times to(c) and ri(c) being defined in (jS.lip and Lemma [331 respectively. 
Then, 

w{T,y) > Wc{T,y) := _^ J,^)^^^,^) Vc,i{t + to{c) - Ti{c),y{l + t)) (4.9) 
for (r,y) € [ti(c),oo) x R. 

Step 2. An explicit family of supersolutions. Let us introduce the following family 
of functions: 

^«(^'y) = (^J (t^-I^IJ ' (4.10) 



We easily obtain by direct calculation that Fji is a classical solution of ([4.3p for y ^ 0, and 
for all parameter values R>0. However, near ?/ = 0, it is only a supersolution both in the 
weak and the viscosity sense. The latter is straightforward to verify using the definition of 
viscosity subsolutions and supersolutions with jets, as in the classical survey [10]. Let us 
mention at this point that these functions can be used in a comparison argument to give 
an alternative proof of (|4.7p . 

Remark 4.1. This family of functions arises naturally if we think about asymptotics. 
Indeed, as already mentioned, we formally expect that the asymptotic profiles of (j4.3p 
should be given by solutions of the stationary Hamilton-, Jacobi equation 

\Vw\P~^ - w = 0, (4.11) 
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supported in some ball B{0,R), that is 

Making then the "ansatz" that, for large times, the solution of ()4.3p should behave in a 
similar way as its limit, we write 

Integrating the resulting ordinary differential equation for C{t), we arrive at the family of 
explicit exact profiles Fn given by (|4.10p . 

Step 3. Constructing suitable subsolutions. We now face the problem of finding 
suitable subsolutions with similar behaviour. Since the -Fr's are classical solutions to (|4.3|) 
except at y = 0, we expect to be able to construct also a family of subsolutions based on 
them. To this end, we consider the "damped" family F^^^^p defined by 

FR^r, y) ■■= ^ \^—^) [ r+l ~ J ^ > y) e Q , (4.12) 

for parameters i? € (0, 1), -i? € (0, 1], and /? € (1/2, 1]. Observe that, since (p— l)/(p — 2) > 
1, Fr^^^p and |VF^,^,^|P-2VFfi,^,^ both belong to Ci([0, oo) x (M \ {0})). For i9 G (0, 1), 
/3 G (1/2, 1], T > and y / 0, we calculate 



drFR,^,p - (ApF^j,^,^ + y ■ VFr^^^p - ^^^^fl,^,/?) + I^Fr^^M^'^ - Fr. 

> ^ (,aP-2 I3{r + R) \ ^l/(p-l) 

(l + ^)2^R,i,/3 l + rV r + 1 



< T?(i - 



L,/3 



1 p-2f(3{T + R) 



1 + T P-IV •^ + 1 



Analyzing the sign of the last expression and taking into account that G (0, 1), we obtain 
that has the following properties: 

Fr,-&,i3 is a classical subsolution to (|4.3|) in ^ATi\ 

{{r, y)eQ : T > r2(i?, /3) , < |y| < A'^,^(r)} ^^'^''^ 

with 
and 

Fr,^,I3 vanishes for |y| > ^ and r > 0. (4-15) 

Let us notice here that both the edge of the support of and the constant Kr^p{t), 

where the behaviour changes, do not depend on {}. While the two properties (|4.13p and 
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()4.15p are suitable for our purpose, the function F^^^^p does not behave in a suitable way 
near y = (where it is a viscosity supersolution) and in an asymptotically small region 
near the edge of its support (where it is a classical supersolution). However, we already 
have a positive bound from below for w in a small neighbourhood of y = by ()4.6p which 
allows us to remedy to the first bad property of i*)?,^,/?. More precisely, we infer from (|4.6|) 
that 

I /3^2\ (p-i)/(p-2) 
w{T,y)>C^:= — (-f\ >0, (T,y) G [0,cxd) x5(0,n/2), 

whence 

w{T,y)>^>FR,^,p{T,y), (T,y) e [0,oo) X S(0,r,/2), (4.16) 

provided that 

< I? < min{l,C5}. (4.17) 

Consider next 

T>T2{R,P) and Kr^p{t) < \y\ < ^^^^f ^ ■ 

Then 



FR,^,^{T,y) < § 



p_ 2\ (P-l)/(P-2) I \ (p-l)/{p-2) 



p-lj \p-2 l + T 



(1 + -r)(p-l)/(p-2) ■ 

Now, if c € (/3, 1), we have 

|y|(l + r) < /3(t + R) < Zc,i,(i_c)/(i+c) + c(r + to(c) - ri(c)) 

as soon as 



(4.18) 



r > r3(c, R, P) := + ^^"^^^^ " " . (4.19) 



In that case. 



r \ 1 r , r \ c ^ n , ^^ -^c,(i-c)/(i+c) 

= il + r)iP-mP-2) ^cAr + ro{c)-n{c),yil + T)) = ^ 

according to the properties (j3.12p of Vc^i- Recalling ()4.9p and (j4.18p we realize that 

i^ii,^,/3(r, y) < Wc{t, y) < w{t, y) , Kr,p{t) < \y\ < ^iZ±3. , (4.20) 

T + 1 

provided 

cE(/3,l), < min{l,A4^(i„,)/(i+c)}, r > max {ri(c), r2(ii, /3), r3(c, i?, /3)} . (4.21) 

After this preparation, we are in a position to establish a positive lower bound for w on 
the bah B{0, 1 — e) for any e € (0, 1/4). Indeed, we fix e G (0, 1/4), choose c = 1 — e, 
i? = /3 = 1 — 2e, and define 

r4(e) := max | ^^^i^^, r2(l - 2e, 1 - 2e),r3(l -e,l - 2e, 1 - 2e) 1 . 
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As Ti{e) > ri(l — e)/e, ()4.8p guarantees that B{0, 1 — 2e) C Vwir^ie)) and there is thus 
rris G (0, 1) such that 

w{T4{e),y)>m,, yeB{0,l-2e). (4.22) 
Now, for S (0, 1) satisfying 

< ^9 < min {m„ C5, Mi_,,,/(2_e)} (4.23) 

we infer from (gUD, ([TO]) . ([CTTll . ([TO]) . ([OT]l . and K22\i that 

Fi_2£,i»,i-2£(r,y) < w{T,y), \y\ € |y , Ki„2£,i-2e(T)| , r > r4(e), 

and 

Fi_2e,^,l-2e{uie),y) < ^ < rUe < w{u{e),y) , y < |y| < i^l-2£,l-2e(T4(e)) < 1 - 2e. 
It then follows from (|4.3|) . (|4.13|) . and the comparison principle that 

T 

Fl-2£,^,l-2e{T,y) < w{T,y) , y < |yl < Ki-2e,l-2e{T) , T > T^ie) . 

Recalhng ()4.15p . (|4.16p . and ()4.20p . we have thus estabhshed that 

Fi-2eAi-2e{T, y) < w{t, y) , T € [Ti{e), oo) X M , (4.24) 
for all I? e (0, 1) satisfying g^S]). 

Step 4. Positive bound from belov^^. For e € (0, 1/4), fix e (0, 1) satisfying ()4.23p . 
According to (|4.24p . we have, for r > r4(e) + 1 and y G 5(0, 1 — 3e), 

, , , /p_2V^"'^/^^"'^ /(l-2e)(r + l-2e) , \ (p-i)/(p-2) 



p- ly V •^ + 1 

- '[p-lj [ r + 1 j + 

„ /2(p-2)e2\(?'-^)/(P-2) 

We have thus proved that, for all e G (0, 1/4), there are fJ-e > and t^{£) := T4(e) + 1 such 
that 

0<^i,<w{T,y), (T,y) G [r5(e), 00) X 5(0,1 -3e). (4.25) 

Step 5. Convergence. Viscosity argument. To complete the proof, we use an 
argument relying on the theory of viscosity solutions in a similar way as in the paper [16] 
for the subcritical case of p.ip with q G {l,p — 1). We thus employ the technique of 
half-relaxed limits [5] in the same fashion as in [19] Section 3] and |16j . To this end, we 
pass to the logarithmic time and introduce the new variable s := log(l + r) along with the 
new unknown function 

w{t, y) = u{\og (1 + r), y) , (r, y) G [0, 00) x M . 
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Then, drw{T,y) = e ^dsUj{s,y) and it follows from ()4.3p and ()4.4p that u solves 

e"^ (^dsU} - ApUJ -yVoj + + |Vtj|P"^ - ^ = , (s, y) G Q , (4.26) 

with initial condition = uq. We readily infer from Lemma |4. II that 

||6j(s)||i + ||l^(s)||oo + ||Vw(s)||oo <C4, s>0, (4.27) 
uj{s,y) = for s > and |y| > 1 + Rie~'' . (4.28) 

We next introduce the half-relaxed limits 

^*{y) ■= liminf uj{X + a,z) and uj*{y):= limsup uj{\ + a,z), 

ia,z,\)^{(T,y,co) (o-,2,A)-i.(o-,j/,oo) 

for (s, y) G Q, which are well-defined according to the uniform bounds in (|4.27|) and indeed 
do not depend on s > 0. Then, the definition of tj* and uj* clearly ensures that 

< uj^iy) < uj*{y) for y G M , (4.29) 

while the uniform bounds ()4.27p and the Rademacher theorem warrant that and lo* 
both belong to W^'°°{R). Finally, by Proposition [O applied to (|4.26p . and oj* are 
viscosity supersolution and subsolution, respectively, to the Hamilton-Jacobi equation 

H{C,VC):=\VC\P-^ -C = in M. (4.30) 

Our aim is now to show that uj^, > tj* in M (which implies that lo^, = uj* by ()4.29p ). Since 
UJ* and are subsolution and supersolution to (|4.30p . respectively, such an inequality 
would follow from a comparison principle which cannot be applied yet without further 
information on uj* and w*. We actually need to prove the following two facts: 

(a) uj^{y) = uj*{y) = if |y| > 1, 

(b) ^^*(y) >L^,(y) >0 if yG 5(0,1), 

and then to follow the technique used in |16| to conclude that = uj* and identify the 
limit. 

To prove assertion (a), let us take y G M with \y\ > 1. We then deduce from ()4.28p 
that there exists si(y) > such that uj{s,y) = for s > si{y). Pick sequences (o"„)n>i, 
(An)n>i, and {zn)n>i such that (Tn 0, Xn ^ oo, Zn y, and uj{an + A„,z„) uj*{y). 
On the one hand, there exists ni[y) > such that (T„ + A„ > si(y) for any n > ni{y); 
hence uj{an + A^, y) = for any n > ni[y). On the other hand, we can write: 

|w((T„ + \n,Zn) - w((T„ + A„,y)| < \y - Z„j||Vw((T„ + A„)||oo < C4\y - Zn\ 0, 

hence uj*{y) = = uj^:{y) for any y G M with \y\ > 1. In addition, since uj* and are 
continuous, it follows that uj* = uj^^ = also for \y\ = 1, hence assertion (a). 
To prove assertion (b), let us take y G i?(0, 1). Then, there exists e G (0,1/4) such 
that y G -8(0, 1 - 4e). Since 1 - 3e > 1 - 4e, there is r2(y) > such that B{y,r2{y)) C 
i?(0, 1 — 3e) and we deduce from (|4.25p that there exists S2(e) := log {T^{e) + 1) > such 
that uj{s,z) > fi^ for any s > S2{£) and z G B{y,r2{y))- We now pick sequences {an)n>i, 
(A„)n>i and {zn)n>i such that 0, A„ oo, Zn y, and w((T„ + A„,2;„) uj^{y). 
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Then there exists again n2{y) > such that (T„ + A„ > S2{y) and z„ € B{y,r2{y)) for any 
n > n2{y)- Consequently w((T„ + A^, z„) > for any n > n2{y)- This readily implies that 

> f^*(y) > ^£ > 0, hence (b) is proved. 
We follow the lines of |T6] and introduce 



WM = W*{y) = ^.;*(y)(P-2)/(p-l)^ (4.31) 

p — 2 p — 2 

for any y G B{0,1). From Proposition I7.2| it follows that W^, and VK* are respectively 
viscosity supersolution and subsolution of the eikonal equation 

|VC| = 1 in 5(0,1), 

with boundary conditions W*{y) = W^{y) = for \y\ = 1 and are both positive in B{0, 1). 
Using the comparison principle of Ishii [13], we find that W*{y) < hence they 

should be equal by (|4.29p . It follows that = u* = W in i?(0, 1), where W is the 
viscosity solution to ()1.13p 

\VW\P-^-W = in B(0,1), W = on dB{0,l), 

which is actually explicit and given by 



(p-l)/(p-2) 



as stated in Theorem ll.il In addition, the equality = uj* and (|4.28|) entail the conver- 
gence of uj{s) as s — )• DO towards W in L°°(M) by Lemma 4.1 in or Lemma V.1.9 in [^. 
We end the proof by rephrasing the two scaling steps and arriving in this way to p.lUp . 

□ 



4.3 Proof of Theorem [IHt N >2 

We now prove Theorem 11.11 for radially non-increasing initial data to the problem posed 
in dimension N > 2. We follow the same steps as in dimension = 1, and we only 
indicate below the main differences that appear. These differences are mainly given by the 
appearance of the new term 

^^^^\drw\P~'^drW, r=\y\, (4.32) 
r 

in the radial form of the p-Laplacian term. As we shall see, performing carefully the same 
steps as for dimension = 1, we find that this term does not change anything in an 
essential way. We follow the same division into steps as the case A^ = 1. 

Step 1. Thanks to the construction performed in Section [3.2| this step is the same as in 
dimension A^ = 1. 

Step 2. Due to the appearance of the extra term ()4.32p in the radial form of the equation 
()4.3p . we check by direct calculation that, in dimension N > 2, the function Fji given by 
formula (|4.10p is now a strict supersolution to (|4.3p in Q. Indeed, for y ^ 0, 

drFR - f ApFfl + y . VF^ - ^^Fn] + \VFr\p-^ -Fr = , , Fr . 

l + rV p-2 j (1 + T)\y\ 
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Moreover, its singularity at y = is now stronger. This seems to introduce a new difficulty, 
but we will see that it can be handled by the same perturbation techniques. Let us notice 
at this moment that can be used for upper bounds in the same way as in the case 
= 1, and that still solves the limit Hamilton- Jacobi equation ()4.1ip . 

Step 3. In order to construct subsolutions starting from the family of functions F/j, we 
follow again the ideas of the case = 1. The calculations will be different in some points. 
We consider again the damped family -Fr,,?,/3 defined in (|4.12p for ii € (0, 1), i9 € (0, 1), 
and (3 G (1/2, 1]. For y 7^ we have 



Y := drFR. 



1 



^F, 



1 + T 

1 + r 



p — 1 



+ 



p-2 
(l+r)|2/| ^'^^P 



Fr. 



■d,l3 



m-1 



Fr. 



^,13 



At this point, we further assume that \y\ > r*/2, the radius being defined in Corol- 
lary [27T1 and that 

(4.33) 



-2< (1-/3K 



2{N - 1) 



Since < 1, we obtain 



Y < ^F],%-'^ 



1 + r 

< ^(1 - , 
from which we conclude that 



:i 



(1 +T)r* 
1 p-2 f pjT + R) 

T+1 



-2^ p(p-2)/(p-l) 



1+r p-1 



is a classical subsolution to (14.31) in 



{(r,2/) € g : r > r2(i?,/3) , (r,/2) < \y\ < Kr^^{t)} 



(4.34) 



where T2{R, f3) and Kr^p{t) are still given by ()4.14p . We now proceed as in the one 
dimensional case to establish (j4.24p for all € (0, 1) satisfying (j4.23p along with 



< 



iV- 1 ' 
for ()4.33p to be satisfied. 

Steps 4 & 5. The final steps of the proof are similar to the one dimensional case. □ 



5 Arbitrary initial data 

So far, we have proved Theorems ll.ll and ll.2l for radially non-increasing initial data satis- 
fying ()1.3p . We now extend these two results to general initial data satisfying ()1.3p . 



Proof of Theorems 11.11 and II. 2L Since uq ^ 0, there are xq € K^, tq > 0, and ryo > 

such that uq{x) > 2rjQ for x € B{xQ,rQ). Then, there exists a radially non-increasing 
initial condition uq satisfying p.3p but with support in B{0,ro) and such that ||uo||oo < 
and uo{x) < uo{x — xq) for x £ M.^ . Similarly, there is a radially non-increasing initial 
condition Uq satisfying (|1.3p but with support in B(0,Rq) for some Rq > Rq and such 



23 



that Uoix) > II 

"^^ollcxD for X G B(0,Rq). Denoting the solutions to (jl.ip by u and U with 
initial conditions uq and Uq, respectively, the comparison principle and the translational 
invariance of (jl.ip ensure that 



u(t, X + xq) < u(t, x) <JJ {t, x) , (t, x) € Q . 

Moreover, since 

(p-2)|x + xo|^(^"'^/^^-') 



(5.1) 



1 



log* 



1 



{p-2)\x\ 
logt 



(p_l)/(p_2) 



< 



(P - l)|3;o| 

logt 



and Theorems 11.11 and 11.21 applv to both u and U, the expected results follow from ()5.ip . 
□ 



Appendix. Some results about viscosity solutions 

We state, for the sake of completeness, some standard results in the theory of viscosity 
solutions, that we use in the proof of Theorem 1 1.101 The first one concerns the "viscosity" 
limit of a family of small perturbations and can be found in [Gj Theorem 4.1]. 

Proposition 7.1. Let be a viscosity subsolution (resp. a viscosity super solution) of the 
equation 

Hs{x,Ue,Vue,D'^Ue) = in , 

where is uniformly hounded in all variables and degenerate elliptic. Suppose that {n^} 
is a uniformly bounded family of functions. Then 

u*{x) := limsup ^^(y) (7.2) 

is a subsolution of the equation 

H4x,u,Vu,D^u) = 0, (7.3) 

In the same way, 

u^(x) := liminf Usiv) 

is a supersolution of H* (x,u, Vti, D^u) = 0. Here, and H* are constructed in the same 
way as and u* . 

In other words, this result can be applied to asymptotically small perturbations of a 
known equation, as we do in Section HI 
We also use the following result: 

Proposition 7.2. Let u € C(il) be a viscosity solution of 

H{x,u,Vu)=0 inn, (7.4) 

where Q C and H is a continuous function. If ^ E C (M) is an increasing function, 
then V = $(u) is a viscosity solution of 

H {x, ^-^{v{x)), {^-^)'{v{x))Vv{x)) = 0. (7.5) 
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The same result holds true for subsolutions and supersolutions and can be found in [6]. In 
particular, we use this result in order to pass from the Hamilton-Jacobi equation | Vti|^~^ — 
u = to the standard eikonal equation |Vu| = 1. Finally, we also use the (now standard) 
comparison principle for viscosity subsolutions and supersolutions of the eikonal equation, 
that can be found in |13j . 
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